Oscillatory Instabilities of Standing Waves in One-Dimensional Nonlinear Lattices 
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In one- dimensional anharmonic lattices, we construct nonlinear standing waves (SWs) reducing to 
harmonic SWs at small amplitude. For SWs with spatial periodicity incommensurate with the lattice 
period, a transition by breaking of analyticity versus wave amplitude is observed. As a consequence 
of the discreteness, oscillatory linear instabilities, persisting for arbitrarily small amplitude in infinite 
lattices, appear for all wave numbers Q 7^ 0, n. Incommensurate analytic SWs with \Q\ > n/2 may 
however appear as 'quasi-stable', as their instability growth rate is of higher order. 

PACS numbers: 63.20.Ry, 45.05,+x, 05.45.-a, 42.65. Sf 



A wellknown and rather spectacular phenomenon oc- 
curing in many nonlinear media (e.g. fluids or optical 
waveguides) is the modulational (Benjamin-Feir) insta- 
bility (MI), by which a travelling plane wave breaks up 
into a train of solitary waves (see e.g. 0]). It is also well- 
known that wave propagation in many continuous non- 
linear media is well described by Nonlinear Schrodinger 
type equations, where the solitary wave trains are de- 
scribed by spatially periodic and stable standing wave 
(SW) solutions, the so-called cnoidal envelope waves . 

An analogous instability of propagating waves occurs 
also in discrete systems (e.g. anharmonic lattices), where 
in the case of soft (hard) anharmonicity, waves with small 
(large) wave numbers are unstable and typically 

break up into arrays of intrinsically localized modes or 
discrete breathers (see e.g. |^||)- However, concerning 
the possible existence and stability of spatially periodic 
SWs in anharmonic lattices, much less has been known 
up to now. Although stable SW solutions of the form 
Acos(Qn) cos(uit) generally exist in harmonic lattices as 
linear superpositions of counterpropagating waves with 
equal amplitudes A and frequencies w but opposite wave 
numbers ±Q, there is a priori no guarantee that these 
SWs will remain stable in the presence of anharmonicity. 

Here, we consider the general class of one-dimensional 
(ID) anharmonic lattices described by a discrete non- 
linear Klein-Gordon (KG) or Schrodinger (DNLS) equa- 
tion. We first propose a method to generate (numeri- 
cally) exact time-periodic SW solutions to these equa- 
tions. These SWs, whose spatial periodicity can be ei- 
ther commensurate (Q/2n rational) or incommensurate 
(Q/27T irrational) with the lattice period, are found as 
'multibreather' solutions Q from the so called 'anticon- 
tinuous' limit of zero inter-site coupling, and they reduce 
to the harmonic SWs for small amplitudes. Then, we in- 
vestigate the stability properties of these SWs, and find 
as a most striking result that for infinite lattices, SWs 
with small but nonzero amplitude are unstable through 
an oscillatory instability for all < Q < ir M. Thus, 
the SWs are unstable also for wave numbers where MI 
for propagating waves does not occur. It is important to 



note that the additional SW instabilities appear as a di- 
rect consequence of the lattice discreteness, and have not 
been found within continuum approximations j§] . More- 
over, the dynamics resulting from the SW instabilities is 
fundamentally different from that of the MI, as the lat- 
ter generally is non-oscillatory (characterized by purely 
imaginary eigenfrequencies of the linearized equations (^) 
below) while the former is oscillatory (corresponding to 
complex eigenfrequencies with nonzero real part) ||. In 
this Letter, our main results will be stated and motivated 
briefly; full technical details will be given elsewhere |l0[ . 

We thus consider a ID chain of classical anharmonic 
oscillators with a general on-site potential V(u) and har- 
monic inter-site coupling Ck > 0, yielding a discrete non- 
linear KG equation for the particle displacements «„, 

u n + V'(u n ) - C K {u n +i + - 2u„) = 0. (1) 

Here, we shall be mainly concerned with the small- 
amplitude dynamics of KG chains with small inter-site 
coupling Ck , which can be well described by a DNLS ap- 
proximation [p|,^| Jl0| (its range of validity is investigated 
in more details in |10|]). For small oscillations, the on-site 
potential V(u) can be expanded as V(u) = u 2 /2 + cra 3 /3 
+ /3u 4 /4 +..., where the linear oscillator frequency is set 
to 1. Linearizing Eq. ([!]) yields propagating or SW solu- 
tions with wave number Q and frequency ujq(Q) given by 
the dispersion relation cjq(Q) = \+ACk sin 2 Q/2. For the 
nonlinear Eq. (|l|), we search for small-amplitude time- 
periodic solutions as u n (t) = J2 P e lpulbt , for which 
the dynamics will be almost harmonic with frequency u)b 
close to u>o(Q) for some Q. Then, J 1 — 1 = 2(5 will be 
of order Ck, and allowing for a slow time dependence 
of the Fourier coefficients offi yields p0[ , at order Ck, 
a DNLS equation for the dominating coefficient aih 1 de- 
scribing the leading-order nonlinear effects. Defining A = 
— 5a 2 /3 + 3/3/2 (so that for soft [hard] potentials, A < 
[A > 0]) and ip n = ^/[Aja^, the DNLS equation reads 

iip„ = Sip n - a\ip n \ 2 ip n + C(ip n+ i + xp n -i - 2tp n ), (2) 

with er = sign(A) and C = Ck/2. We will assume a = — 1 
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without loss of generality, since changing ip n — > (— l)™i/>*, 

<5 — > 4C — (5 is equivalent to changing the sign of er. 

(v) 

For time-periodic solutions all a n T , and therefore ?/> n , 
will be time-independent and can be chosen real. Then, 
vanishing the left-hand side of Eq. (||) defines a nonlinear 
symplectic map S : (ip n ,ip n -i) —> {ipn+i,i>n) m tne 2D 
real plane, which, with the change of scale ip n — > VCip n , 
depends only on the parameter 5' = 8/C. As is well- 
known [ jTl| ], this map exhibits a rich variety of orbits, 
including elliptic and hyperbolic fixpoints and periodic 
cycles, KAM tori, Aubry-Mather Cantor sets (Cantori) 
p2]| and chaotic orbits. Searching for nonlinear SWs in 
continuation of the linear SWs ip n = Acos(Qn + <p>), 
with 5 = 4Csin 2 Q/2 = Sq(Q), only the periodic and 
quasiperiodic orbits which can be continued to zero am- 
plitude by varying 5' are of interest. For small ampli- 
tudes, these orbits must be located close to the elliptic 
fixpoint F = (0,0), and Q must be close to its linear 
rotation number 6(6') given by 6 = 5o{0)- As a < 0, the 
rotation number for orbits rotating around Fq increases 
with increasing radius, and since the map S as well as its 
orbits at fixed rotation number depends continuously on 
6' , SWs with wave number Q exist only for 6 < 5q(Q). 

When Q = 2nr/s (r and s irreducible integers), there 
are two families of commensurate SWs, each of them 
represented by one of the two s-periodic cycles (here- 
after called 'h- cycle' resp. 'e- cycle) which bifurcate in 
pair from F at 6 = S (Q) ||-||. The h-cycle is hy- 
perbolic for all 8' , while the e-cycle is elliptic for S close 
to 5q(Q) and hyperbolic with reflection below a critical 
value S = S C (Q). When Q/2ir is a generic irrational num- 
ber, the trajectories representing incommensurate SWs 
emerge from Fq as the KAM torus with rotation number 
Q. At the critical value 5 — 6 C (Q), the KAM torus breaks 
up in a transition by breaking of analyticity |l^] and bi- 
furcates into a Cantorus and its associated 'midgap' tra- 
jectory |14|. Thus, for S < S C (Q) the Cantorus and the 
midgap trajectory define two families of non-analytic in- 
commensurate SWs, which merge into the unique ana- 
lytic SW denned by the KAM torus for S > 5 C (Q). 

To give a global representation of these SWs, we con- 
sider the large-amplitude limit <5 — > — oo, or, equivalently, 
the anticontinuous limit C — for fixed S < 0. From this 
limit, where stationary solutions to Eq. (||) can only take 
the values ip n = ±v^5 and ip n = 0, we can describe the 
SWs as multibreathers 13 5| identified by a coding se- 
quence a' n defined as ip„ = y/—6o~' n . By continuing the 
SW map orbits to large negative 5' , we find their cod- 
ing sequences as a' n = Xo(Q n + 4>), where Xo(x) is the 
27r-periodic odd function defined for x € [0, tt] as 



Xo(a;) 



1 for (n-Q)/2<x< (tt + Q)/2 
elsewhere 



(3) 



Q/2w). We will call these SWs 'type H\ and they have 
the property (of importance for the stability analysis be- 
low) that their coding sequences do not contain any con- 
secutive +1 or —1. For the particular phases <j) — <j) m , 
x = Qm + <j) is at a discontinuity point of Xoi x )) so that 
a' n has two consecutive +1 (or —1). These SWs, called 
type E, correspond to e-cycles or to midgap trajectories. 

The representation (^) is useful also for numerically 
calculating directly the exact SWs of the original KG 
equation (|l|), using standard methods (e.g. Newton 
schemes) for continuation of multibreathers from the an- 
ticontinuous limit. Continuing these solutions versus u>b 
to zero amplitude, we typically find (Fig. [I]) a transition 
by breaking of analyticity of the hull function xi x ): de- 
fined as u n (0) = x(Q n + 0), for incommensurate SWs in 
the KG chains as well as in DNLS chains. 




For all phases <f> ^ 4> m = ±(ir — Q)/2 — mQ (m inte- 
ger), this describes SWs corresponding either to the h- 
cycles (for rational Q/2ir) or to the Cantori (for irrational 



FIG. 1. Hull functions for SWs in a KG chain (gj) with 
V(u) = | {e~ u - if, C K = 0.05, Q/2-k = 21/55 ~ (3— y/E)/2 
= (7 q and (a) u>b = 1.072 resp. (b) ujb = 1.060. The breaking 
of analyticity occurs at ujb — 1.068 for Q/2n — gq. 

We now turn to the investigation of the dynamical 
stability properties of these SWs. As before, we con- 
sider small-amplitude solutions to the KG equation (|l|) 
with small Ck, so that the DNLS approximation (||) is 
well justified. With the substitution ip n — > ip n + e n (t), 
where ip n is real and time-independent and e n is small, 
the linearization of Eq. (||) yields a standard Hill equa- 
tion for e„. Its eigenmodes can be found by setting 
e»(t) = \{u n + v^e^ + i« - <)e- iu ' f , where the 
eigenfrequencies u>i are given by the eigenvalue problem: 

C V n = (2C - S ~ 1pl)v„ - C(v n+ i + Vn-i) = LUlU n 

C\u n = (2C - S - 3ipl)u n ~ C(u n+ i + Un-i) = uJiv n . (4) 

The SW ip n is linearly stable if and only if all eigenfre- 
quencies uji are real. 

At the anticontinuous limit (8 < 0, C = 0), the eigen- 
values of (^) will be located either at uji — (corre- 
sponding to sites with codes \a' n \ — 1) or at wj = ±5 
(corresponding to a' n = 0). For commensurate SWs 
(Q = 2-Kr/s) the lattice periodicity of ip n in (H) implies 
that the degenerate eigenvalues will form bands as C is 
increased from zero (see Fig. ^). There will be s—2r pairs 
of bands originating from u>i — ±5, and r pairs of bands 
originating from uji — (which always remains in the 
spectrum). For incommensurate SWs, the band struc- 
ture will be Cantor-like. In all cases, bands originating 
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from ±S can be shown to move initially only along the 
real axis when C is increased from zero, while the be- 
haviour of the bands originating from will be different 
for E and H type SWs. As the coding sequences for E 
type SWs contain pairs of consecutive +1 (or —1) and the 
anharmonicity is soft, it follows from a general result || 
that some eigenvalues always move out on the imaginary 
axis as soon as C ^ 0, making these SWs dynamically un- 
stable for 5 < and small C. Moreover, we find JhJ that 
this instability persists for all S < Sq{Q), so that SWs of 
type E are unstable for any nonvanishing amplitude. By 
contrary, for H type SWs the theory of effective action 
H can be used to show that all eigenvalues initially will 
move along the real axis, and thus the stability of these 
SWs is preserved for C not too large (Fig. || (a)). 

Increasing 8 (for fixed C) the bands broaden, and the 
main gap separating the two classes of bands originat- 
ing from resp. ±<5 shrinks. At some value 5 = Sk{Q) it 
vanishes, and, as shown in Fig. |^(b), eigenvalues from the 
two classes of bands will collide and move out in the com- 
plex plane (since they have opposite Krein signature ^|). 
Thus, oscillatory instabilities appear for the SW. These 
instabilities have also been confirmed by direct numerical 
Floquet analysis of the KG Eq. (0) for small C K ©• 

0.2 
0.1 

lm((0|) 
-0.1 
-0.2 





(a) 







-10 

Fig. 2(a), Morgarile, Phys. Rev. 
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Fig. 2(b), Morgarile, Phys. Rev. Lett. 

FIG. 2. Eigenvalues of Eq. tyty for a type H SW with 
Q = 3tt/4 and (a) 8' = -8.0 resp. (b) 5' = -7.5 (C = 1). 



For 6 > Sk(Q) bands continue to overlap, generating 
new instabilities. In many cases, the SW remains unsta- 
ble until it vanishes at 8 = 5q(Q), although for some Q 
with small s the stability is temporarily regained for some 
interval of S where no bands overlap. But for all commen- 
surate SWs, there is a final interval 5i(Q) < 5 < 5q(Q) 
where bands with opposite Krein signature overlap, and 
as a consequence all commensurate SWs are unstable for 
small but nonzero amplitude. This can be proven | p0| by 
considering the linear limit (6 = 5o(Q),ipn = 0)> where 
the spectrum of (^) is given by lui = ±2C\vQ(q)\ with 
vq(o) = cosQ — cosg. In this limit, eigenvalues with 
opposite Krein signature (given by sign(i/g(g)) [^0)) will 
overlap in an interval around ll>i — 0. Since for com- 
mensurate SWs a finite number of bands with nonzero 
widths appear for 5 < 6q(Q) when ^ 0, a continuity 



argument |T(| yields that bands with opposite Krein sig- 
nature will continue to overlap close to lji = and cause 
instabilities until 5 reaches some value Sx(Q) < 5q(Q). 
(In the linear limit, the SW is stable in spite of the band 
overlap, since the eigenmodes are uncoupled.) 

By contrast, for small-amplitude (analytic) incommen- 
surate SWs the Cantor-like nature of the spectrum of (Q) 
(with 'bands' of zero width) prevents the use of a similar 
argument for proving instability. However, the analysis 
outlined below shows that also these SWs are unsta- 
ble, but the instabilities may be of higher order and thus 
more difficult to detect numerically or experimentally. 

The instability of analytic SWs is proven using the 
method of 'band analysis' ||). Its basic idea is to embed 
the eigenvalue problem (|J) for a non-Hermitian operator 
into a wider eigenvalue problem for a Hcrmitian operator, 
which we choose as (with u> a real parameter) 



-UJ 



— LJ 

Co 



{U n } 
{V n } 



= E 



{u„} 



(5) 



Then, each eigenvalue E v is a smooth, real function 
('band') of u>. The real eigenfrequencies uoi of (Q) are 
given by the intersections of these bands with E = 0, 
and an instability occurs whenever a band loses a pair of 
intersections with E = 0, i.e., when a 'gap' opens around 
E = 0. 

In the linear limit (ijj n — > 0, S — > 5q(Q)) the bands can 
be labelled by a wave number q and describe straight 
lines E±(q;uj) — 2G(cosQ — cosq) ±uo. Thus, for each q 
there are two bands (doubly degenerate since E±(q; lo) = 
E±(—q;uj)) 1 and bands with different q intersect when 
E^(q-;u) = E + (q + ;u). When tp n describes an analytic 
SW with small amplitude A, tp^ can be treated as a per- 
turbation and expanded as 4>n = ^2 P f p e lp ( 2nQ+ ^ , where 
the coefficients f p are of order A 2 ^ for \p\ > 1, and thus 
exponentially decaying for p — > oo. This introduces a 
coupling between eigenvectors of (^) with wave numbers 
q and q + 2pQ with coupling strength proportional to f pi 
and gaps of width AE cx 2\f p \ will open at each inter- 
section where E_(q;uj) = E + (q + 2pQ;uj). If these gaps 
surround the axis E = 0, instabilities occur for the SW. 

For fixed p these intersections describe an ellipse in the 
(lu, i?)-plane given by 



2C sinpQ 



E p - 2CcosQ 
2C cospQ 



= 1. 



(6) 



Numerical calculations of the set of bands E{uj) confirm 
the existence of these ellipses of gaps (see Fig. ||) . The el- 
lipse with p = 1 corresponds to the largest gaps (and thus 
the strongest instabilities if the gaps surround E — 0) and 
is for 5 = 5 (Q) tangent to E — at uj — 0, lying below 
E = for 7r/2 < Q < n and above for < Q < n/2. 
Considering also the shift of the bands caused by the 
static term /o, we find |10| that when < Q < it/ 2 gaps 
will open around E = in the lower part of the ellipse, 
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and thus these waves become immediately unstable to first 
order (as for commensurate SWs with Q ^ it ]Tc|]). 



E 




Fig. 3, Morgante, Phys. Rev. Lett. 0) .— . 

FIG. 3. Band spectrum of (g) for a SW with 
Q/2% = 34/89 ~ era {Q > tt/2) and 5' = 3.4 (C = 1). The 
visible gap openings occur along the ellipse ([]) with p = 1. 

By contrast, for irrational Q with 7r/2 < Q < 7r all 
first-order gaps open strictly below E = (Fig. ||), so 
that no first order instabilities develop for these SWs. 
However, as is seen from (^), the ellipse corresponding to 
the pth order intersections will intersect E = if | cospQ| 
> | cosQ|, and therefore there is always ap such that the 
SW is unstable to order p for arbitrarily small amplitude. 
But as the width of the gaps and the maximum instabil- 
ity growth rates are proportional to \f p \, the instabilities 
become very weak when the smallest p yielding an insta- 
bility becomes large (for Q close to tt). We might there- 
fore view these analytic SWs as 'quasi-stable'. A typical 
example of the dynamics resulting from this instability in 
a KG chain is shown in Fig. [|, where, after the initial os- 
cillatory dynamics and an intermediate regime of mainly 
translational motion, an apparently chaotic final state 
appears. Further details concerning the long-time dy- 
namics of the unstable SWs will be published elsewhere. 




Fig. 4, Morgante, Phys. Rev. Lett. 

FIG. 4. Time evolution of an analytic SW (at integer multi- 
ples of the original period) , perturbed only by numerical trun- 
cation errors, in a KG chain ([l]) with Ck = 0.03, uit = 1.044, 
and other parameters as in Fig. [I]. 

In conclusion, we have shown that for a general class of 
anharmonic lattices, different types of SWs with spatial 
period commensurate or incommensurate with the lat- 



tice exist, but are generically unstable through oscillatory 
instabilities for small amplitudes. Thus, in discrete sys- 
tems, the nonlinearity-induced coupling between counter- 
propagating waves leads to wave break-down for arbitrar- 
ily weak anharmonicity even if the individual propagating 
waves are modulationally stable. The general nature of 
this result suggests that these instabilities should be ob- 
servable in macroscopic as well as microscopic contexts. 
In particular, as was recently experimentally verified [ 51| , 
nonlinear optical waveguide arrays provide a direct appli- 
cation of the DNLS equation, and thus these systems are 
also good candidates for detection of SW instabilities. 
M.J. acknowledges a EC Marie Curie fellowship. 
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